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1 INTRODUCTION 

In quantum theory, exact statements on solutions of eigenvalue equations — describing 
bound states — may be obtained by some rather simple and elementary methods. Here, 
these methods are reviewed and illustrated at the simplest nontrivial (semi-)relativistic 
wave equation, the spinless Salpeter equation, which represents a well-defined standard 
approximation to the famous Bethe-Salpeter formalism. For more details, see Ref. 

1.1 Relativistic Wave Equation: Spinless Salpeter Equation 

1.1.1 Bound States in Quantum Field Theory: Bethe— Salpeter Formalism 

The appropriate framework for the description of bound states within some relativistic 
quantum field theory is the Bethe-Salpeter formalism. There a bound state |M) (with 
momentum K and energy E) of, say, fermion and antifermion (with masses M\, M 2 
and momenta Pi, P 2 , resp.) is represented by its Bethe-Salpeter amplitude defined, 
in momentum space, in terms of the time-ordered product of the field operators of the 
two bound-state constituents between vacuum and bound state by 

*(P) = exp(iKX CM ) J d 4 X exp(iPX) (0|T(V>i(Xi) $ 2 (X 2 ))\M(K)) , 

with the total momentum K = Pi + P 2 , relative momentum P, center-of-momentum 
coordinate Xcm, and relative coordinate X = X\—X 2 of the bound-state constituents. 
The Bethe-Salpeter (BS) amplitude ^ satisfies the Bethe-Salpeter equation 

S^(P 1 )^(P)S 2 \-P 2 ) = 7^4/ d 4 Q/C(P,Q)*(Q) , 

which involves two dynamical ingredients: 

• The full fermion propagator Si(P) of some particle i (i = 1, 2) of momentum P 
and mass Mj is usually approximated by its free form S~q(P) = — i (7^ P^ — Mi), 
with Mi interpreted as some effective ("constituent") mass and the propagator 
understood as an effective one. 

• The BS kernel K(P, Q) is defined perturbatively as sum of all two-particle (BS-) 
irreducible Feynman diagrams for two-particle into two-particle scattering (and 
thus gauge-dependent!); it is given, e.g., in lowest order QED from one-photon 
exchange (in Feynman gauge) by 

The BS equation is formally exact. Nevertheless, it faces some well-known problems: 

• There is no means to compute the BS kernel beyond perturbation theory. 

• In general, it is not possible to find the exact solutions of the BS equation (except 
for very few special cases, like the famous Wick-Cutkosky model for two scalar 
particles interacting by exchange of some massless scalar particle). 

• In non-Abelian gauge theories like QCD, free propagators are incompatible with 
a confining BS kernel because Dyson-Schwinger equations connect propagators 
and kernel. 
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The usual way out is the reduction of the BS equation by a series of approximations: 
1. Eliminate any dependence on timelike variables: 

• Adhere to the static approximation to the BS kernel, i.e., assume that the 
kernel depends only on the relative t/iree-momenta of initial and final state: 
/C(P, Q) = /C(P, Q). (This is tantamount to ignoring all retardation effects 
by assuming instantaneous interactions.) 

• Define the equal-time wave function $(P) = J dP ^(P, Pq). 
This leads to the Salpeter equation 

A+7o/C(P,Q)$(Q) 7 oA2 



r d u C 
J (27 



d 3 g 

(2tt) 3 



Ar7o/C(P,Q)$(Q)7oA 2 f 



Ml 



P + y/p? + Mj 2 + y^l + Mf 
with the positive/negative-energy projection operators for particle i (i — 1,2) 



AT = 



_ v /p|+M?±7o(7-P, + M,) 



2^P 2 + M 2 



The Salpeter equation is the equation of motion for the equal-time wave function 
$ with full relativistic kinematics but in static approximation for the kernel K. 

2. Neglect the second term on the r.h.s. of the Salpeter equation by assuming 

E - y/Pl + Ml - yfpl + M| < E + y/Pl + Ml + yJP% + Ml . 

This leads to the reduced Salpeter equation 
[e - VP? + Ml - yjj* + M|) $(P) = / ^ A+ 70 /C(P, Q) $(Q) 7 o A 2 " . 

3. Neglect all the spin degrees of freedom of the involved bound-state constituents. 

4. Assume that the BS kernel /C(P, Q) depends only on the difference of the relative 
momenta, which means that it is of convolution type: /C(P, Q) = /C(P — Q). 

5. Restrict the whole formalism exclusively to positive-energy solutions ip. 

Under these simplifying assumptions and approximations, the BS equation reduces to 
the spinless Salpeter equation 



P? + Ml + y/pi + Mf + V(X) 



ip = Etp , 



involving an interaction potential, V(X), arising as the Fourier transform of /C(P — Q). 
In the center-of-momentum frame of the two bound-state constituents, i.e., for K = 0, 
this equation reads H\tp) = E\ip), with the Hamiltonian 



H = JP 2 + Ml + JP 2 + Ml + V(X 
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1.1.2 Equal-Mass Case 

For equal masses of the two bound-state constituents (i.e., assuming Mi = M 2 = M), 
the two-particle spinless Salpeter equation is equivalent to its one-particle form. To see 
this, consider the two-particle Hamiltonian H with an interaction represented, e.g., by 
a central potential of power-law form: 

H = 2 VP 2 + M 2 + ^ k n R n , i2=|X|. 

One may always perform a scale transformation of the phase-space variables X, P by 
some arbitrary scale factor A, 

p = AP , x = — , 

since this rescaling preserves the fundamental commutation relations: [x, p] = [X, P]. 
Fix A to A = 2, which implies 

H = V4P 2 + 4M 2 + ]T k n R n = ^p 2 + 4M 2 + £ k n 2 n r n , 

neZ neZ 

identify the one- and two-particle mass and coupling-strength parameters according to 
m = 2 M and a n = 2™ k n , n G Z, and arrive at the equivalent one-particle Hamiltonian 

H = \jv 2 + ™ 2 + J2 a n r n , r = l x l ■ 



1.1.3 Prototype: (One-Particle) Spinless Salpeter Equation 

In view of the above, it is sufficient to consider the self-adjoint Hamiltonian H = T+V, 
with the "square-root" operator of the relativistic kinetic energy of a particle of mass m 
and momentum p, 

T = T(p) = vV + m 2 , 

and arbitrary coordinate-dependent, static interaction-potential operators V = V(x). 
The spinless Salpeter equation is then nothing else but the eigenvalue equation for 
H, i.e., H\xk) = E k \xk), k — 0, 1, 2, . . ., for Hilbert-space eigenvectors \xk) and energy 
eigenvalues 

E ^ (Xk\H\xk) 
(Xk\Xk) 

It thus represents the simplest relativistic generalization of the Schrodinger equation of 
standard nonrelativistic quantum theory. N.B. The semirelativistic Hamiltonian H is 
a nonlocal operator, i.e., 

• either T in configuration space 

• or, in general, V in momentum space 

is nonlocal. Therefore, it is difficult to obtain rigorous analytic statements on solutions 
of this equation of motion. 
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1.2 Example: The (Spinless) Relativistic Coulomb Problem 

For illustrative purposes, let us consider a (spherically symmetric) Coulomb potential, 
with interaction strength parametrized by some (dimensionless) coupling constant a: 

V(x) = V c (r) = -— , r = Ixl , a > . 
r 

For this case, some important pieces of knowledge have been accumulated until now: 
• An examination of the spectral properties of H with a Coulomb potential reveals 



the essential self-adjointness of if (i.e., its closure is self-adjoint) for a < — ; 

2 

the existence of the Friedrichs extension of if up to a critical value a c = — ; 

7T 

a strict lower bound on the ground-state energy E (that is, on if) , given by 



E > m d 1 — I I for a < — . 

In this sense, the Coulombic Hamiltonian if is a reasonable operator up to a c . 
For part of the allowed range of a, an improved lower bound may be derived H]: 



1 + Vl-4q2 1 

tin> m\ for a < - . 

- V 2 2 

Lower bounds on the spectrum of H may be found numerically with the help of 
the generalized "local-energy" theorem ||: Assume that 

1. the Fourier transform V(p) of V(x) is strictly negative, except at infinity, 
as is the case for an attractive Coulomb potential, 

2. the spectrum of H is discrete, and 

3. the ground state of H exists. 

Define the "local energy" 

J d 3 gy(p-q)0(q) 
with some suitably chosen, positive trial function 0(p) > 0. ThenF] 



5(p) = T(p) + 



inf£(p) < E Q < sup£(p) . 
p P 

This theorem restricts the ground-state energy Eq{o) as a function of a to some 
remarkably narrow band. In particular, at a c one finds 

E 

0.4825 < — < 0.4842910 for a = a c . 
m 

Consequently, E (a = a c ) is definitely nonvanishing. 
However, even for the Coulomb potential the eigenvalues of H are not known exactly! 
1 The lower bound even holds if the spectrum is not purely discrete. 
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2 Analytic Upper Bounds on Energy Levels Q5j 

We are interested in exact analytic expressions for upper bounds on the eigenvalues of 
H. This only makes sense if the operator under consideration is bounded from below. 
Consequently, assume that the potential V is such that H is bounded from below. For 
the Coulomb potential, this has been rigorously proven 0. 



2.1 Bounds on Eigenvalues: Minimum— Maximum Principle 

The primary tool for the derivation of rigorous upper bounds on the eigenvalues of a 
self-adjoint operator is the well-known minimum-maximum principle. Among several 
equivalent formulations the most convenient one for practical purposes is the following: 
Let if be a self-adjoint operator bounded from below with eigenvalues Ek, k = 0, 1, . . ., 
ordered according to E < E\ < E 2 < . . ., and let D d be some d- dimensional subspace 
of the domain of H . Then the fcth eigenvalue E k (when counting multiplicity) satisfies 

E k < sup 7., r 7 for fc = 0,l,2,... . 

2.2 Operator Inequalities 

We would like to replace the problematic kinetic-energy square-root operator in H by 
some more tractable operator. One way to achieve this is to use the min-max principle 
in order to compare eigenvalues of operators: 

• Assume the validity of a (generic) operator inequality of the form H < O. Then, 

^ mom 

E k < SUp — 7T7TT — < SU P 



• Assume that D k +i is spanned by the first k + 1 eigenvectors of O (corresponding 
to the first k+1 eigenvalues E , Ei, . . . , E k of O if the latter are ordered according 
to E < %_ < E 2 < . . .). Then 

sup ... = E k . 

|^>GD fe+1 WW) 

Every eigenvalue of H is bounded from above by the related eigenvalue of O: E k < E k . 



2.3 The "Schrodinger" Bound 

The simplest upper bound on the eigenvalues of H may be easily found by exploiting 
the positivity of the square of the (self-adjoint, since T is self-adjoint) operator T — m: 

0< (T - m) 2 = T 2 + m 2 - 2mT = p 2 + 2m 2 ~ 2mT . 

Assuming m to be positive, this may be converted into an operator inequality for T, 
T < m + p 2 /(2m), which entails an operator inequality for H: H < m + p 2 /(2 m) + V. 
For the Coulomb potential, the required Schrodinger energy eigenvalues are given by 

with the total quantum number n, expressed in terms of the radial quantum number n T 
(n r = 0, 1, 2, . . .) and the orbital angular momentum £ (I — 0, 1, 2, . . .) by n — 
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2.4 A Straightforward Generalization 

Real improvement is made by considering the positivity of the square of the (obviously 
self-adjoint) operator T—fi, where /i is an arbitrary real parameter (with the dimension 
of mass): 



implies a set of operator inequalities for T (see also Ref. 0), T < (p 2 + m 2 + / u 2 )/(2 //) 

V n > 0. This translates into an operator inequality for H: H < (p 2 +m 2 +/i 2 ) / (2 fj)+V 

V fi > 0. Hence, according to the min-max principle, E^ is bounded from above by the 
eigenvalue Es^(n) of the Schrodinger-like Hamiltonian on the r.h.s. of this inequality, 
Ek < Es J k{/J') V fi > 0, and thus also by the minimum of these Schrodinger-like bounds: 



For the Coulomb potential, the corresponding "Schrodinger" energy eigenvalues are 



with the total quantum number n = n T +£+ 1. Minimizing i?s,n(/-0 w.r.t. fi then yields 



These bounds hold for all a < a c and arbitrary levels of excitation, and they improve 
definitely the Schrodinger bounds, for any n: 



(For fi = m one necessarily recovers the Schrodinger case.) The comparison of these 
analytic (!) upper bounds with their numerically obtained counterparts shows that the 
relative error of these analytic bounds (for a < 0.5) is for the ground state (n r = £ = 0) 
less than 4.5 % and for the level n r = 0, £ = 1 less than 0.1 %. 

2.5 Rayleigh— Ritz Variational Technique 

An immediate consequence of the min— max principle is the famous Rayleigh 
Ritz technique: Restrict the operator H to the subspace Dd by orthogonal projection P 
onto Dd'. H := PHP. Let Ek, k = 0,1, ... , d—1, denote all d eigenvalues of H, ordered 
according to E < E\ < . . . < E d -\. The kth. eigenvalue (counting multiplicity) of H, 
Ek, then satisfies Ek < Ek, k = 0, 1, . . . , d — 1. 

Now, if Dd is spanned by linearly independent basis vectors \ipk), k — 0,1, ... , d—1, 
the eigenvalues E are simply determined by diagonalizing the dxd matrix ((tfji\H\ipj)), 
i,j — 0,1, ... ,d — 1, i.e., as the d roots of the characteristic equation 



0< (T - fi) 2 = T 2 + fj, 2 -2fiT = p 2 + m 2 + fj 2 -2fiT 



E k < mm E S)k (fi) . 





min£s,„Gu) < E s , n for a ^ . 



det((^|#|^>-£(^;>) = , i,j = 0,1,... ,d-l 
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2.6 Variational Bound for the Ground State 

As a first test, let us apply the Rayleigh-Ritz variational technique to the ground state, 
by considering the case k — (i.e., an only one- dimensional subspace). In this case, the 
min-max principle reduces to 



E < 



(</#> 



(i.e., the ground-state energy E is less than or equal to any expectation value of H). 
However, one may consider simultaneously more than one one-dimensional trial spaces 
in order to obtain an optimized upper bound according to the following prescription: 

1. Choose some suitable set of trial states {|?/>a)} (with elements distinguished from 
each other by some variational parameter A). 

2. Calculate all the expectation values E(\) = (tp\\H\ip\) / (ipx\ipx) of H w.r.t. 

3. Determine (from the first derivative of E{\) w.r.t. A) that value A min of A which 
minimizes E(X). 

4. Compute E(X min ) (i.e., the minimal expectation value of H in the Hilbert-space 
subsector of the chosen trial states). 

The outcome of this procedure is an optimized upper bound on Eq, viz., Eq < E(X m i n ). 

In order to get rid of the troublesome square-root operator in H, we adopt a trivial 
but fundamental inequality for the expectation values of a self-adjoint operator O and 
its square w.r.t. any state in the domain of O: 



o\ ., < 



ie> 2 



\i m) ' 

We consider this inequality, of course, for T, and apply the resulting inequality to H: 

^\T + V\^) 



E < ^^^2ZL < 



<V|T 2 |V> + (1>\V\1>) - 



<?/#) -\ m) (m \ <v#> m) 



+ m l + 



For the Coulomb potential, a special choice for the coordinate-space representation 
of trial vectors is suggestive: the well-known (normalized) hydrogen-like trial functions 

^a(x) = W — exp(-Ar) , A>0. 



This leads to a one-parameter set of upper bounds, E < a/A 2 + m 2 — a A for all A > 0, 
with absolute minimum E < m a/1 — « 2 , which is identical to the previous generalized 
(operator) bound for n = 1 and, for a ^ 0, lower and thus better than the Schrodinger 
bound 

£s,o = m ^1 - y^j . 

Thus, the variational technique yields indeed improved upper bounds on energy levels. 
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2.7 Energy Levels at the Critical Coupling Constant of the 
Relativistic Coulomb Problem 

In order to find, for the Coulomb potential, upper bounds on the energy levels at a c , 
we employ basis vectors \ipk) (labelled by a positive integer k = 0, 1, 2, . . .) of Dd, given 
in configuration-space representation by 



(2 m) 



2fc+2/3+l 



\ AnT(2k + 2p + 1) 



r k+ ^ exp( 



-mr 



(3 > , m > 



The parameter (3 allows, for a given a, of the complete cancellation of the divergent 
contributions to the expectation values of the kinetic energy T for large momenta p and 
the Coulomb potential Vq(t) at small distances r: (3 = (3(a). (3 is implicitly determined 
as a function of a, e.g., for the ground state by a = (3 cot(/3 7r/2). Consequently, a c is 
approached for (3 — > 0. For our \4>k)i singularities arise only in the ground-state matrix 
elements, i.e., in (ip \T\ip ) and ( , 0o|Vc(r')|V'o)- It is a simple task to evaluate all relevant 
matrix elements for arbitrary (3. 

For simplicity, introduce a dimensionless energy eigenvalue e by E =: (2/7r)me. 
The resulting characteristic equation for (3 = is typically of the form 





1 41n2 


-2-s 


v/2 
3 


£ 


_ \ 




det 


V2 


e 


17 


- e 




= 




3 


"71 


15 ~ 







V 



The roots of this equation (which may be calculated algebraically up to d = 4) read, for 
d = 1, e = 2 (2 In 2 — 1), entailing the ground-state upper bound E /m = 0.4918 . . ., for 

d = 2, 

e = — (60 In 2 - 23 ± ^(60 In 2) 2 - 4800 In 2 + 1649 

entailing the ground-state upper bound E /m = 0.484288 . . ., while, for d = 4, some 
lengthy expression implies the ground-state upper bound E Q /m = 0.4842564 . . .. From 
d = 2 on, all these analytic bounds lie well within the numerical "local-energy range." 



2.8 Generalized Laguerre Basis 

Upper bounds on eigenvalues may be improved by enlarging to higher dimensions d 
or by spanning it by a more sophisticated set of basis states. A rather popular choice of 
trial functions involves the generalized Laguerre polynomials L^\x) for parameter 7: 
these are orthogonal polynomials defined by the power series 

tf'M = g(-i>'(^)£ 

and normalized, with the weight function x 7 exp(— 2), according to 

OO 

|dxx^exp(-x)4 7) W4' ) (^) - £ + S *k> ■ 
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To construct a trial vector corresponding to a state with orbital angular momentum £ 
and its projection m, introduce two variational parameters: fi (with dimension of mass) 
and f3 (dimensionless). This state is defined by the configuration-space representation 



(2/i) 2 ^ +1 fc! i+n-i , x r (2e+2(3) (? ]v (Q 
\}T(2£ + 2{3 + k + l) r 6XP( ^ T)Lk ( 2 ^)y^> 



with the orthonormalized spherical harmonics ye m (fl), depending on the solid angle Q. 
Normalizability of the trial states requires \i > and (3 > —1/2; ipk,im{^) then satisfies 
the (standard) orthonormalization condition / d 3 a; ^ <m (x) ipk'j'm'fe) = $kk' $te' $mm> ■ 

2.8.1 Power-Law Potentials 

Let us investigate interaction potentials of power-law form: V(r) = J2n a nr bn , r = |x|, 
where the real constants a n and b n are only constrained by the semiboundedness of H: 
b n > — 1 if a n < 0. The matrix elements of V are easily worked out: 



(A\V\^j) = V r (2£ + 2 3 + i + l)T(2£ + 2 3+ i + 1) ? (2 u) K S ^ 



r(2* + 2/3 + i + l)r(2* + 2/3 + j + l) V (2//) 6 ™ ^ 0s f tls\ 

r(2£ + 2/? + 6 n + t + s + 1) . 



i + 2£ + 2/?Wj + 2£ + 2/3 
* - t ) \ j - a 



2.8.2 Analytically Evaluable Special Cases 
Orbital Excitations 

Consider just orbital excitations by restricting to i = j = but allowing for arbitrary I. 
Only for definiteness, fix 3 to 3 = 1. In this case, the matrix elements of T are given by 

with the hypergeometric series 

„, N T(w) ~ r(u + n)r(w + n) z n 

-F(m, f ; u>; 2) := ^ r . 

T(u) T(v) ^ T(w + n) n\ 

There are several possibilities to get rid of F: In the ultrarelativistic limit (m = 0), e.g., 
the kinetic-energy matrix element becomes 

For instance, for the linear potential V(r) = ar, a > 0, minimizing (ip \H\ifj ) w.r.t. 
\x leads to the minimal upper bound 

r(£ + 2) 



min^ol-f^lV'o) = 2 v2 a 



At>0 



r(£ + |) • 



In the limit £ — > 00, this reduces to 



l i m(min^ | Wo )) 2 = 8a(£+^ 

describing linear Regge trajectories, i.e., [_E(£)] 2 cx £, in striking accordance with other 
findings based on different considerations |7|, || . 
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Radial Excitations 

Consider pure radial excitations by focusing to states with I = 0. The matrix elements 
of T then typically take the form 



r(2/3 + i+l)r(2/3 + j + l) 2tt 
i + 2(3 \ f j + 2(3 ' 



^ f ("2) 4+s 

t=0 s=0 



t!s! 



3 — s 



T(f3 + t + l)T((3 + s + l)L, 



ts ) 



with J ts , for /i = m, and for 2 /3 integer (thus non-negative), i.e., 2 /? = 0, 1, 2, . 
by 



'2p + t + s+\t-s\ + l 



-1 |t-s| 

E 

n=0 



\t - S\ 

n 



X !( ^J r 2 008 It 



T[2f3 + t + s + - 



-1 2/3+t+s+2 

E 

n=0 



2/5 + t + s + 2 
n 



x [M'^^ , ]r(2/? + t + s + i 



2 J cos It - 



given 



By this, analytic expressions for the matrix elements (ipi\H\ipj) of H with a power-law 
potential may be found. Up to d = 4, these matrix elements are, at least in principle, 
algebraically accessible. For d > 4, the energy matrix [{ip i \H\'ipj)) must be diagonalized 
numerically, without, however, the need of (time-consuming!) integration procedures. 
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